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PREFACE 


This  study  has  as  its  goal  the  development  of 
procedure  for  correcting  measured  Hall  carrier  profiles, 
which  have  been  distorted  from  the  true  carrier  profiles  by 


a  surface  depletion  effect.  All  the  necessary  lab  work  was 
carried  out  at  the  Avionics  Laboratory. 

1  wish  to  express  great  appreciation  for  the 
assistance  and  guidance  of  many  individuals  who  assisted  me 
during  this  study.  I  give  special  thanks  to  Or.  Y.  K. 

Yeo,  my  thesis  advisor,  for  his  invaluable  guidance.  I 
also  give  heart  felt  thanks  to  C.  Geesner,  who  deposited 
the  encapsulants,  and  J.  Ehret  for  performing  the  ion 
implantation. 
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y  ABSTRACT 

V  method  was  developed  for  correcting  a  carrier 
profile  by  Hall  Measurements  for  the  surface  depletion 
effect.  The  method  assumes  that  the  experimental  data  is 
to  be  fit  with  a  Pearson  type-IV  curve. 

In  order  to  correct  the  measured  profile  it  is 
necessary  to  know  the  depletion  width  after  each  etch.  The 
depletion  widths  are  dependent  upon  the  real  profile  and 
not  the  measured!  therefore  if  one  only  has  measured  data, 
the  depletion  width  cannot  be  directly  determined.  In 
order  to  gain  insight  into  the  problem  a  procedure  called 
"reverse  correction"  is  developed.  It  creates  a  measured 
profile  from  a  real  profile. 

The  method  for  correction  is  the  following.  First 
assume  a  real  profile  is  known.  Calculate  the  depletion 
widths  from  this  profile.  From  these  depletion  widths  and 
the  Pearson  type-IV  fit  to  the  original  measured  data 
construct  a  new  assumed  real  profile.  Repeat  this 
procedure  self-consistently  until  convergence  of  the 
assumed  real  data  to  a  stable  profile.  This  is  the  true 
profile.  Convergence  is  expected  in  approximately  7 
iterations. 


SURFACE  DEPLETION  CORRECTION  TO  CARRIER  PROFILES 


BY  HALL  MEASUREMENT 


I.  INTRODUCTION 


In  recent  years,  the  use  of  ion  implantation  as  a 
means  of  introducing  impurities  into  a  semi -insulating 
semiconductor,  for  the  purpose  of  device  fabrication,  has 
been  given  considerable  attention.  For  instance,  active 
layers  in  microwave  FETs  (Field  Effect  Transistors)  must  be 
precisely  controlled  to  depths  of  about  1  micron  (1:14). 
This  requirement  is  difficult  to  meet  by  standard 
techniques  of  epitaxial  growth  or  thermal  diffusion. 
However,  ion  implantation  allows  one  to  control  the  doping 
profile  by  simply  varying  the  energy  and  fluence  of  the 
incident  ions.  Thus,  theoretically,  a  continuous 
variation  of  the  fluence  and  energy  can  produce  any 
distribution  desired  (2:9).  In  addition,  ion  implantation 
has  been  found  to  produce  good  uniformity  of  implants  over 
large  regions,  as  well  as  good  reproducibility  from  wafer 
to  wafer  (3:122). 

Of  particular  interest  is  silicon-implanted  GaAs, 
since  Si  is  a  relatively  light  atom  and  is  therefore 
readily  accelerated,  with  modest  energy  requirements,  to 
velocities  necessary  to  produce  useful  implant  depths. 
Furthermore,  it  has  found  that  low  dose  implants  of  Si  in 


r 


GaASf  have  produced  high  electrical  activation 
efficiencies,  which  is  more  difficult  to  achieve  with  other 
n-type  dopants  such  as  S  and  Se  (4:27). 


There  are  many  other  applications  for  ion  implantation 
which  show  considerable  promise.  For  instance,  in  a 
research  project  conducted  for  NASA,  Woo  (5)  constructed 
complementary-metal -oxide-semi conductor  / 

silicon-on -sapphire  circuits  (CMOS/SOS)  by  ion  implantation 
which  were  one-third  faster  than  those  fabricated  by 
diffusion.  Finally,  other  researchers  have  demonstrated 
the  feasability  of  producing  electrically  insulating  layers 
in  device  materials  by  ion  in^lantation,  in  lieu  of  the 
more  common  mesa  techniques  (1:17). 

In  general,  for  device  applications,  it  is  necessary 
to  accurately  measure  the  carrier  profiles  of  implanted 
semiconductors .  The  two  common  electrical  means  for 
accomplishing  this  are  capacitance -voltage  (C-V)  profiling 
and  secondly  by  the  differential  Hall  method.  The  former 
makes  use  of  a  reversed  biased  Shottky  barrier  to  produce  a 
variable  depletion  width,  determined  by  measuring  its 
capacitance.  As  the  reverse  bias  is  increased  the 
depletion  width  extends  deeper  into  the  material, 
decreasing  the  capacitance,  and  uncovering  more  dopant 
ions.  The  density  of  carriers  at  the  edge  of  the  depletion 
width  is  thus  related  to  the  differential  change  of 
capacitance  with  bias  voltage.  There  are  two  main 


disadvantages  associated  with  C-V  profiling.  The  first  is 
that  the  profile  cannot  be  measured  inside  the  initial 
depletion  width*  which  occurs  at  zero  bias  voltage.  The 
second  major  problem  is  that  C-V  profiling  is  useful  only 
in  n-type  materials. 

Differential  Hall  profiling,  on  the  other  hand,  may  be 
used  in  both  n-  and  p-type  materials.  This  technique  is 
characterized  by  measuring  of  the  sheet  carrier 
concentration  in  carriers/cm  ,  which  corresponds  to  the 
total  number  of  electrically  active  carriers  below  the 
surface  per  unit  surface  area.  Next,  a  thin  layer  of  the 
implanted  material  is  removed  by  a  chemical  etch.  A 
subsequent  measurement  of  the  sheet  carrier  concentration 
reveals  the  density  of  activated  ions  in  the  layer  which 
has  been  etched.  An  obvious  disadvantage  with  this 
technique  is  that  it  is  destructive,  and  therefore, 
measurements  may  not  be  repeated.  Another  serious 
disadvantage  is  that  a  surface  depletion  region,  associated 
with  the  filling  of  acceptor-type  surface  electronic  states 
(SES),  results  in  the  measurement  of  a  profile  which  is  not 
only  shifted  from  the  true  profile  but  also  contracted. 

This  translation  is  due  to  the  fact  that  Hall  measurements 
make  use  of  only  mobile  carriers,  thus  only  the  carriers 
past  the  edge  of  the  depletion  width  are  measured. 
Furthermore,  the  distortion  occurs  because  the  depletion 
widths  are  a  function  of  carrier  concentration,  which,  of 


course,  varies  with  depth.  Small  depletion  widths  occur 
where  the  concentration  of  implanted  ions  is  large  and 
large  depletion  widths  where  the  concentration  is  small. 

THESIS  OBJECTIVE  AND  SCOPE 

It  is  the  objective  of  this  research  to  develop  a 
method  for  correcting  the  Hall  measured  profile,  thereby 
obtaining  a  profile  which  accurately  represents  the 
distribution  of  free  carriers.  The  corrected  profiles  will 
then  be  compared  to  the  profiles  obtained  by  C-V  profiling. 

The  substrate  material  used  was  <100>  oriented  undoped 
semi -insulating  GaAs.  Ion  implantation  was  carried  out  at 
room  temperature  at  an  ion  energy  of  100  keV.  The 
substrates  were  implanted  with  doses  of  2,  4,  6,  and 
8  X  10^^  as  well  as  1  and  3  x  10^^  cm”^ .  The  last 
dose  produces  degenerate  conditions  in  the  sample.  The 
methodology  developed  for  correcting  for  correcting 
non-degenerate  doses  must  be  easily  adapted  to  the  case  of 
degenerate  doses. 

SEQUENCE  OP  PRESENTATION 

Chapter  II  provides  background  information  on 
ion-implantation,  with  an  emphasis  on  ion  range  statistics. 
In  addition  this  chapter  explains  the  theoretical  basis  of 


surface  electronic  states,  and  establishes  a  relationship 
for  the  depletion  width.  Chapter  III  presents  the 
experimental  procedures  and  reviews  the  methods  of  Hall  and 
C-V  profiling.  With  this  preparation,  chapter  IV  details 
the  reasons  for  the  distortion  of  the  Hall  profile  as  well 
as  the  method  for  correcting  it.  A  theoretical  profile  is 
used  to  demonstrate  the  effect  of  surface  depletion  on  the 
measured  profile.  Finally,  Chapter  V  presents  the 
experimental  results,  as  well  as  conclusions  concerning  its 


success . 


II.  BACKGROUND 


RANGE  DISTRIBUTIONS 


The  penetration  depth  o£  an  incident  ion  in  a  solid 
material  is  governed  by  various  factors,  among  which  are 
the  mass  and  energy  of  the  ion,  the  mass  of  the  target 
atom,  the  temperature  of  the  target  material,  etc.  The 
range  can  be  accurately  determined  if  one  possesses  a 
thorough  knowledge  of  the  stopping  power  (-dE/dx)  or 
specific  energy  loss.  This  differential  function  is 
dependent  on  the  instantaneous  energy  of  the  ion  and  has 
two  main  contributions,  the  electronic  and  nuclear  stopping 
powers,  that  is. 


dE/dx»  (dE/dx)  +  (dE/dx)  (1) 

e  n 

Electronic  stopping  is  characterized  by  small  angular 
elastic  collisions  between  the  ion  and  the  target  atom 
electrons  resulting  in  the  ionization  of  the  material.  The 
nuclear  stopping  term  is  due  to  the  collisions  of  the  ion 
with  the  nucleus  of  the  target  atom,  where  the  field  is 
non-coulombic  and  "Rutherford  type"  scattering  takes  place. 
Nuclear  energy  loss  dominates  at  low  ion  energies  while 
electronic  loss  is  prevelant  at  large  ion  energies 


Lindhard,  Sharff,  and  Shiott  (LSS)  have  made  numerical 
calculation  of  the  ion  range  in  amorphous  materials  by 
using  the  Thomas-Fermi  interatomic  potential  to  calculate 
the  contribution  from  the  nuclear  energy  loss  term.  In 
addition,  they  have  computed  a  mean  square  fluctuation  in 
the  energy  loss  for  each  collision,  resulting  in  a  mean 
square  fluctuation  in  the  total  range  ^  .  This 
calculation  of  straggling. assumes,  as  a  first  approximatior 
a  Gaussian  range  distribution  (2:36-39).  That  is,  the 
distribution  of  implanted  ions  with  depth  is  given  by 


N(x) 


-(x-R  )2 


Where  R  is  the  projected,  or  mean  range  perpendicular  to 

tr 

the  surface,  ^  is  the  fluence  or  dose  in  ions /cm  ,  and 
<7  is  the  standard  deviation  about  the  projected  range. 

For  instance,  for  an  energy  of  100  keV,  the  project  range 
of  Si  in  GaAs  is  0.085  itm,  with  a  standard  deviation  of 
0.0442  ftm.  Figure  1  illustrates  the  LSS  distribution. 

However,  experimental  evidence  has  shown  that 
depending  upon  the  relative  importance  of  the  nuclear 
collision  process,  the  distribution  can  be  quite  asymetric, 
necessitating  higher  order  moments  of  the  distribution. 
Another  factor  influencing  this  asymmetry  is  ion 
channelling,  in  which  incident  ion  travelling  parallel  to 
or  scattered  in  the  direction  parallel  to  a  crystal 


4, 


direction  have  a  projected  range  much  larger  than  expected 
(2:40).  A  distribution  function  which  incorporates  higher 
order  moments  is  the  Pearson  type-IV,  which  satisfies  the 
following  relation 

df  ,  (x-a)f(x) 

h©  +  b^x  +  (3) 

along  with  the  condition 

oo 

yf(x)dx  =  1 

-oo 

The  four  constants  a,  and  b2«  are  functions 

of  the  first  four  moments  of  the  distribution >  that  is,  the 
projected  range,  the  standard  deviation,  skewness,  and 
kurtosis.  According  to  Hofker  (6:46-47),  profiles  obtained 
from  boron  implanted  silicon,  at  various  energies  ranging 
for  30  keV  to  800  keV,  indicate  the  Pearson-IV  distribution 
to  be  far  superior  to  the  Gaussian  distribution  in  fitting 
this  experimental  data. 

Equation  (3)  may  be  integrated  to  obtain  an  analytic 
expression  for  the  distribution.  A  more  detailed  account 
of  the  Pearson-IV  distibution,  including  the  effect  of  the 
projected  range,  standard  deviation,  skewness,  and  kurtosis 
on  the  shape  of  this  distribution  is  given  in  appendix  A. 


SURFACE  ELECTRONIC  STATES 


The  theory  of  rectification  of  a  metal  semiconductor 
contact  developed  by  Shottky  implied  that  the  maginitude 
and  polarity  of  the  contact  potential  would  be  quite 
sensitive  to  the  work  function  of  the  particular  metal  used 
in  forming  the  junction  (7:4).  However,  experimental 
evidence  shows  that  the  contact  potential  is  virtually 
unaffected  by  the  work  function  of  the  metal.  In  1947 
Bardeen  (7:5)  proposed  that  the  potential  barrier  accross 
the  contact  was  produced  by  surface  states,  rather  than  the 
contact  potential.  The  theoretical  existence  of  surface 
states  had  been  surmized  earlier  by  Tamm  and  Shockley, 
however  there  impact  upon  the  electrical  properties  of 
semiconductors  was  not  realized  until  Bardeen's  insight. 

The  existence  of  surface  states,  as  explained  by 
Shockley  (8:317),  is  due  to  the  sudden  asymmetrical 
termination  of  the  periodic  potential  of  the  crystal 
lattice  at  the  surface,  as  shown  in  figure  2.  Thus,  the 
periodic  boundary  conditions  applied  in  Kronig-Penny  theory 
no  longer  hold  at  the  surface.  Shockley  has  shown  that 
this  gives  rise  to  quantum  states,  localized  at  the 
surface,  which  have  energies  within  the  forbidden  gap. 
Furthermore,  it  is  also  appropriate  to  view  the  surface 
electronic  states  as  being  due  to  unfilled  orbitals,  or 
dangling  bonds.  For  instance,  Massies  et  al.  (9:640)  have 


Figure  2.  Asymetrlc  potential  at 
the  crystal  surface 


Figure  3.  Band  bending  at  the  surface 


observed  a  transition  from  a  Ga  3d  to  a  Ga  dangling  bond  on 
the  (100)  face  of  GaAs  (9:640).  Both  Tamm  and  Shockley 
have  shown  there  to  be  approximately  one  surface  state  for 
each  surface  atom. 

The  charging  of  the  surface  states  causes  the  valence 
band  maximum  and  conduction  band  minimum  to  rise  as  they 
approach  the  surface,  as  shown  in  figure  3.  The  depletion 
width  extends  from  the  surface  until  the  depth  where  the 
bands  coincide  with  their  bulk  positions.  As  shown  by 
Chandra  et  al.  (10:646),  the  depletion  width  is  given  by 


In  these  equations,  is  the  difference  in  the 

conduction  band  and  Fermi  level  at  the  surface  in  eV,  and 

is  the  density  of  states  of  states  at  the  conduction 

band  minimum,  and  N,  are  the  concentratitions  of 
a  a 

donor  and  acceptor  dopants,  respectively.  The  value  of  ^ 

B 

for  the  (100)  face  of  n-type  GaAs  has  been  measured  by 
Massies  et  al.  (9:641)  as  well  as  Pianetta  et  al.  (11:617), 


and  has  been  found  to  be  approximately  0.6  eV. 

In  the  case  of  degenerate  statistics,  equation  (5)  is 

no  longer  applicable,  as  the  Fermi  level  moves  above  the 

conduction  band  minimum.  Following  Hattori  (12:6907),  the 

approximation  that  V__  =  when  N(x)  >  N  ,  will  be 

00  o  c 

adopted.  That  is,  in  degenerate  conditions,  the  depletion 
width  extends  approximately  to  the  point  where  th  Fermi 
level  coincides  with  the  conduction  band  minimum. 


ii: 


SAMPLE  PREPARATION 


The  substrate  material  used  in  this  study  was  <100> 
oriented  undoped  semi -insulating  GaAs.  The  2  inch  circular 
wafers  were  cut  into  1/2  inch  by  1/2  inch  square  samples 
which  were  prepared  for  implantation.  This  preparation 
included  cleaning  with  basic-H,  de-ionized  water, 
trichloroethylene,  acetone,  methanol,  and  then  drying  with 
nitrogen  gas.  When  dark  field  microscopy  had  revealed  that 
this  cleaning  had  removed  all  foreign  material  from  the 
substrate  surface,  the  substrate  was  free-etched  with  an 
H2SO^:30%H2O2*H2O  solution  in  a  7:1 »!  ratio  by 
volume  for  3  minutes.  Implantation  of  Si  was  carried  out 

at  room  temperature  with  an  incident  energy  of  100  kev  for 
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doses  ranging  from  2  x  10  through  3  x  10  cm  . 

The  ion  beam  was  directed  7*  off  the  <100>  crystal  axis  in 

order  to  minimize  the  effect  of  ion-channeling.  After 

implantation  the  samples  were  soaked  in  an  HCl  solution  for 

1  minute  in  order  to  remove  a  natural  oxide  layer  which  had 

grown  since  implantation.  The  samples  were  then 

immediately  capped  with  a  1000  A  layer  of  ^^3^4  using 

a  plasma  enhanced  deposition  system.  The  samples  were 


then  annealed  at  850 *C  for  15  minutes  in  flowing  hydrogen 
gas  to  activate  the  implanted  ions.  This  combination 


produces  high  activation  efficiencies#  as  shown  by  Kim 

(4:27),  who  achieved  activativation  efficiencies  of  at 

12  -2 

least  80%  for  6  x  10  cm  implants  of  silicon  in 
GaAs.  After  annealing,  the  encapsulants  were  removed  with 
a  48%  Hydroflouric  acid  solution,  which  required 
approximately  10  minutes.  Each  1/2  inch  by  1/2  sample  was 
then  cut  into  4  quarter  inch  by  quarter  inch  samples. 

Samples  from  each  dose  were  then  profiled  using  the 
capacitance-voltage  method  described  below.  Subsequently, 
several  samples  from  each  dose,  with  satisfactory  C-V 
profiles,  were  chosen  and  sheet  concentration  measurements 
were  made.  Prior  to  sheet  measurements  indium  contacts 
were  placed  on  the  corners  of  the  samples,  and  they  were 
then  annealed  at  375 *C  for  3  minutes  in  order  to  produce 
ohmic  behavior.  Hall  profiles  were  obtained  from  samples 
with  representative  activation  efficiencies  from  each  dose. 


ELECTRICAL  MEASUREMENTS 


CAPACITANCE-VOLTAGE  PROFILING 


A  rectifying  metallic  contact  on  an  n-type 
semiconductor  produces  a  space  charge  region,  or  depletion 
region,  in  the  semiconductor.  Electrons  from  this  region 
are  transferred  to  the  metal,  in  order  that  the  Fermi  level 
be  continuous  accross  the  junction.  If  an  additional 
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reverse  bias  voltage  is  applied  to  the  contact,  the 
depletion  width  may  be  extended  deeper  into  the  material. 
The  space  charge  region  may  be  treated  as  a  voltage 
dependent  parallel  plate  capacitor,  with  capacitance 


where  x  is  the  depletion  width  in  cm,  A  is  the  area  of  the 
2 

contact  in  cm  ,  €  is  the  permittivity  of  the  material  in 

Farads /cm,  and  C  is  the  capacitance  of  the  contact  in 
Farads.  If  the  reverse  bias  voltage  is  increased  by  dV, 
the  space  charge  region  extends  deeper  into  the  material  by 
dx,  which  correspondingly  reduces  the  barrier  capacitance, 
and  uncovers  N(x)Adx  ions  at  the  edge  of  the  depletion 
width.  Henisch  (13:216)  has  related  the  differential 
change  in  capacitance  with  bias  voltage  to  the  carrier 
concentration  at  the  edge  of  the  depletion  width: 

dV 

C-V  profiles  were  obtained  with  the  aid  of  the 
Hewlett-Packard  Model  4061  A  Semiconductor /Component  Test 
System.  The  sample  to  be  profiled  was  placed  implanted 
face  down  onto  two  ducts .  The  ducts  were  then  evacuated 
allowing  mercury  to  rise  to  the  top  of  each,  thereby  coming 


<•> 


into  contact  with  the  implanted  surface.  The  larger  duct 
was  connected  to  the  high  side  of  the  bias  supply  and  the 
small  duct  was  connected  low  side.  Bias  voltages  from  0  V 
to  -8  V  were  then  incrementally  applied  to  measure  the 
profiles.  The  bias  voltage  was  supplied  with  the 
Hewlett-Packard  4140  B  DC  Voltage  source.  The  capacitance 
was  measured  at  each  voltage  step  using  a  0.01  V  1  MHz 
test  signal  supplied  by  the  4275A  Multi -Frequency  LCR  meter 
which  then  measured  the  complex  impedance  and  thereby 
determined  the  capacitance  of  an  assumed  equivalent 
circuit.  This  circuit  was  comprised  of  the  capacitance  of 
the  space  charge  region  in  series  with  the  resistance  of 
the  material  between  the  two  mercury  contacts. 

HALL  MEASUREMENTS 


While  C-V  profiling  makes  use  of  static  charge  to 
determine  the  carrier  profile.  Hall  measurements  use  only 
the  mobile  carriers.  Unfortunately,  the  measurement  of 
Hall  voltages  only  reveals  the  total  number  of  carriers 
associated  with  entire  implanted  region.  Therefore,  in 
order  to  obtain  a  depth  profile  of  carriers,  it  is 
necessary  to  combine  Hall  measurements  with  etching. 

Hall-effect  and  sheet  resistivity  measurements  were 
made  using  the  standard  van  der  Pauw  technique.  This 
technique  requires  that  only  four  contact  be  placed 


anywhere  on  the  periphery  of  a  homogeneously  thick  simply 
connected  sample  of  arbitrary  shape.  However,  in  order  to 
correct  for  this  arbirtary  geometry  L.  J.  van  der  Pauw 
(14:6)  has  shown  that  it  is  necessary  to  measure  the 
voltage  accross  two  adjacent  contacts  when  a  current  is 
passing  through  the  other  two,  as  shown  in  figure  4a,  for 
square  geometry.  In  this  manner,  we  define  the  resistance 
as  Similarly,  is  defined  from  figure 

4b.  From  these  two,  the  sheet  resistivity,  .in  Q/d  ,  is 
given  by. 
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where,  f  is  a  geometrical  correction  factor,  which  is  close 
to  unity  for  ~  1*  Th®  sheet  resistivity  is 

defined  more  accurately  by  the  two  relations,  as  in  Mayer 
(15:4077), 

^  ~  T*  n .  ii  •  d . 

d  1.  (9) 

i 


from  which, 
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where  d.  is  the  thickness  of  the  ith  layer,  /i .  is  the 
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Six  configurations  for  Hall-Effect/ 
sheet  resistivity  measurements 


mobility  in  the  ith  layer,  is  the  carrier 

concentration  in  the  ith  layer  and  d  is  the  total  thickness 
of  the  implanted  region.  Equations  (9)  and  (10)  call 
attention  to  the  fact  that  as  the  carrier  concentration 
varies  with  depth,  so  does  the  mobility,  which  is  only 
natural,  since  the  amount  of  scattering  will  be  a  function 
of  the  impurity  concentration.  In  equation  (9).  parallel 
conductances  at  varying  depths  are  added,  to  obtain  the 
equivalent  conductance  from  the  implanted  region.  Equation 
(9)  is  then  used  to  obtain  an  equivalent  resistance  for  the 
implanted  region. 

The  sheet  Hall  coefficient,  in  ii?cm  /(Vs),  is 
defined  by. 
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where  B  is  the  normal  component  of  the  applied  magnetic 
field  in  Gauss,  and  dV  is  the  change  in  voltage  of 
position  (e)  in  figure  4.  when  the  field  is  applied,  and  I 
is  the  applied  current  in  Amperes.  The  effective  sheet 
concentration,  in  carriers /cm  ,  is  given  by 
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According  to  Petritz  (16:1259).  the  sheet  Hall 


coefficient  is  given  by. 


Petritz  shows  that  this  expression  is  correct  to  the  first 
order  in  B.  A  second  order  term,  which  takes  into  account 
circulating  currents,  is  negligible  for  the  magnetic  field 
used  in  this  study,  0.5  Tesla  or  5,000  Gauss. 

Thus,  by  combining  Hall-effect/sheet  resistivity 
measurements  with  etching  we  can  obtain  the  mobility  and 
carrier  concentration  within  the  etched  interval,  that  is, 

H -  (14 
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The  quantities  (R_)j  and  (  />-)i  are  the  sheet  Hall 
coefficient  and  sheet  resistivity,  respectively,  which  are 
measured  following  the  ith  layer  removal.  Equations  (14) 
and  (15)  are  easily  obtained  from  the  summations  (10)  and 
(13). 

Successive  layers  from  the  the  implanted  region  are  removed 
using  an  etching  solution  of  H2SO^:H202sH20  in 
the  ratio  1:1:100  by  volume  at  0*C.  Measurements  were 
continued  until  all  the  active  layers  were  removed.  Black 
wax,  placed  over  the  contacts  prior  to  the  first  etch,  and 
removed  after  the  last  etch,  revealed  the  total  etched 
depth.  The  etched  depth  was  measured  using  the 
Sloan-Dektak  IIA  Surface  Profiling  System,  and  from  it  an 
etching  rate  was  computed. 

A  block  diagram  of  the  automated 
Hall-ef fect/sheet-resistiviy  measurement  system  is  shown  in 
figure  5.  It  consists  of  a  Digital  Equipment  PDP  11-03 
computer,  four  Keithley  Model  610  electrometers  serving  as 
unit  gain  amplifiers,  a  Keithley  Model  725  programmable 
current  source,  a  Keithley  Model  6900  digital  multimeter 
used  to  measure  voltages,  and  a  Keithley  Model  616  digital 
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Figure  5.  Block  diagram  of  Automated  Hall-Effect/ 
sheet  resistivity  measurement  system. 


electrometer  which  measures  the  current.  As  shown  in 
figure  6,  a  six-position  rotary  switch  is  used  to  measure 
all  the  relevant  parameters  in  each  of  the  configurations 
of  figure  4. 


IV.  THEORY 


THE  PROBLEM 


Equations  (14)  and  (15)  which  give  the  carrier  and 
mobility  profiles  assume  that  the  carriers  which  contribute 
to  the  summations  (9)  and  (13)  extend  from  the  surface 
through  the  end  of  the  implanted  region.  This,  however, 
does  not  take  into  account  the  effect  of  surface  states, 
which  create  immobile  ionized  charges  in  the  depletion 
region.  Therefore,  these  carriers  do  not  contribute  to  the 
Hall  voltages  and  sheet  resistances.  In  fact,  only  the 
mobile  carriers  beneath  the  depletion  width  are  detected. 
That  is,  more  properly,  summations  (9)  and  (13)  should  be 
taken  from  the  edge  of  the  depletion  width,  w^,  through 
the  end  of  the  implanted  region.  After  an  etch,  there  is  a 
new  depletion  width,  w^,  which  extends  from  the  new 
etched  surface,  as  shown  in  figure  7.  Although  the  actual 
measurements  were  made  at  the  depths  of  w^  and  w^,  this 
fact  was  usually  ignored,  and  the  carrier  concentrations 
were  simply  calculated  on  the  basis  of  etched  depth, 
without  taking  into  account  the  surface  depletion  width. 
When  this  surface  depletion  width  is  not  accounted,  the 
measured  profile  obtained  may  be  quite  different  than  the 
true  profile. 

Thus,  we  distinguish  between  the  real  and  the  apparent 


profiles,  the  latter  being  obtained  from  equation  (15). 

From  equation  (14)  it  is  obvious  that  the  apparent  mobility 
profile  does  not  depend  upon  the  etched  depth,  and  hence  it 
is  simply  the  real  mobility  profile  translated  towards  the 
surface  by  the  depletion  width: 


/i^(x  +  w)  » 


(16) 


However,  the  relationship  between  the  real  and 
apparent  carrier  profiles  is  not  as  simple.  In  fact,  the 
d^  in  the  denominator  of  equation  (15)  should  be  replaced 
by  d^  +  (Wj^  "  1^*  Here,  d^  is  the  depth  of  the 

1th  etch,  w^  is  the  depletion  width  after  the  ith  etch, 
and  w^  ^  is  the  depletion  width  after  etch  (i-1) .  This 
correction  acknowledges  the  fact  that  the  depletion  width 
changes  with  each  etch  due  to  the  non-uniform  implant. 
Thus,  the  measured  concentration  in  equation  (15)  is 
related  to  the  real  concentration  by 


n  (x4w) 
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or 
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where  w  is  taken  to  be  the  average  of  the  ith  and  the 
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(i-l)st  depletion  width.  A  differential  relation  may  be 
obtained  by  allowing  d^^  to  approach  zero,  that  is: 

n  (x) 

V  ® 

n  (x  +  w)  *  - 

^  1  +  aw(x)  (18) 

ax 

Equation  (18)  indicates  that  the  true  profile  may  be 
obtained  from  the  measured  profile  if  the  function  w(x)  is 
known.  This  function,  however,  is  determined  by  the  real 
profile,  where  only  the  measured  profile  is  known.  To  gain 
insight  into  the  problem,  it  is  instructive  to  calculate  an 
apparent  profile  from  a  true  one,  or  in  other  words  to 
"reverse  correct". 

REVERSE  CORRECTION 

If  the  real  profile  is  known,  the  measured  profile  may 
be  calculated  with  the  aid  of  (18),  (4),  and  (5). 

Equations  (4)  and  (5),  which  apply  to  a  constant 
concentration,  are  adapted  to  depth  dependent 
concentrations  by  taking  the  average  concentration  within 
the  depletion  width.  Thus,  the  depletion  width  can  be 
self-consistently  determined  in  conjunction  with  equations 
(4)  and  (5).  Specifically,  the  average  carrier 
concentration  within  the  ith  depletion  width  is  given  by 


where  is  the  distance  etched  from  the  surface,  and 
is  an  assumed  depletion  width  from  the  ith  etched 
surface.  The  average  value  is  then  inserted  into  equations 
(4)  and  (5)  to  obtain  a  new  estimation  of  the  depletion 
width  w^.  Equation  (19)  is  re-evaluated  for  the  new 
yielding  an  updated  average  concentration  and 
depletion  width.  This  procedure  is  repeated  until 
convergence  to  the  actual  ith  depletion  width.  This 
usually  requires  approximately  5  iterations  for  an  accuracy 
of  one  thousandths  of  the  depletion  width.  With  all  the 


Wj^'s  calculated  the  measured  profile  may  be  obtained  from 
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where 


Xf  -  (x^  +  x^_j^)/2 


The  reverse  correction  was  carried  out  for  the 
theoretical  LSS  distributions  corresponding  to  the  doses 
2,  4,  6  and  8  xl0^^/cm^  and  1  xl0^^/cm^.  The 


results  are  shown  in  figures  8  through  12.  From  these  we 


CONCENTRATION,  cm' 


make  several  observations  concerning  the  effect  of  surface 
depletion  on  the  measured  profile.  First,  the  peak  of  the 
distribution  is  shifted  in  towards  the  surface,  since  the 
edge  of  the  depletion  width  passes  it  much  earlier  than 
expected.  Secondly,  all  the  measured  distributions  are 
contracted  with  respect  to  the  real  profiles,  for  basically 
the  same  reason.  Third,  the  distortion  for  lower  doses  is 
much  more  significant  than  that  for  higher  doses.  And 
finally,  the  trailing  edge  of  each  measured  profile  drops 
very  sharply.  The  reason  for  this  effect  may  be  seen  in 
figure  7  which  shows  the  etched  surfaces  and  their 
corresponding  depletion  width,  for  a  dose  of 
6  X  10  /cm  .  As  seen  in  the  last  few  etches,  the 
depletion  width  is  increasing  rapidly,  however  the 
concentration  within  adjacent  depletion  widths  is 
decreasing  much  more  rapidly.  This  leads  to  the  sharp 
trailing  edges  shown  in  figures  8  through  12.  Such 
behavior  is  observed  in  profiles  measured  with  the  Hall 
method,  while  C-V  profiles  show  a  longer  tail  in  the 
distribution. 

DEPLETION  WIDTH  CORRECTION 

Reverse  correction  is  a  straight  forward  procedure, 
since  we  have  the  true  profile,  from  which  the  depletion 
widths  are  generated.  Yet,  in  actual  experimental  work,  we 


only  have  the  measured  profile.  The  procedure,  to  be 
developed  below,  for  calculating  the  true  profile  from  the 
measured  one,  shall  be  known  as  "surface  depletion  width 
correction" . 

Equation  (18)  indicates  that  if  the  depletion  width 
after  each  etch  is  known  then  the  problem  is  solved.  A 
first  order  guess  for  these  depletion  widths  may  be 
obtained  by  assuming  a  real  profile,  and  calculating  from 
it  the  depletion  widths  for  the  measured  profile.  Using 
equation  (18),  a  first  approximation  to  the  real  profile 
may  be  obtained.  This  profile  is  certainly  closer  to  the 
true  profile  than  is  the  measured,  since  no  depletion 
widths  were  used  in  the  calculation  of  the  measured 
profile.  This  first  approximation  is,  of  course,  not  the 
true  profile.  It  may  be  checked  for  closeness  by  "reverse 
correction".  If  the  "reverse  corrected"  and  measured 
profiles  coincide  then  the  assumed  real  profile  is,  in 
fact,  the  real  profile.  If  the  measured  and  "reverse 
corrected"  profiles  do  not  sufficently  agree,  then  a  new 
set  of  depletion  widths  are  obtained  from  the  latest 
assumed  real  profile  and  an  updated  real  profile  is  again 
created.  Convergence  of  this  procedure  is  expected  on  the 
basis  that  the  assumption  of  an  incorrect  real  profile  will 
cause  a  perturbing  effect  on  the  next  calculated  profile 
towards  the  true  profile.  The  assumed  real  profiles  will 
continually  change  until  the  only  possible  stable  situation 


IS  obtained,  in  which  one  profile  agrees  with  next 
calculated  profile.  This  means  that  the  calculated 
depletion  widths  are  self-consistent  and  thus  do  not  change 
if  the  reverse  correction  procedure  is  repeated.  A  reverse 
correction  of  the  final  profile  should  prove  this  to  be  the 


case. 


The  correction  of  actual  measured  profiles  was 
obtained  with  the  aid  of  a  computer  program  in  BASIC 
listed  in  appendix  B.  The  procedure  for  correction  is  the 
following.  The  measured  data  is  fit  with  a  Pearson-IV 
curve,  the  parameters  of  which  constitute  the  measured 
profile.  As  a  first  guess  to  the  real  profile,  the 
measured  profile  is  shifted  in  depth  by  the  depletion  width 
corresponding  to  the  maximum  in  its  concentration.  From 
this  assumed  real  profile  an  updated  profile  is  then 
created.  This  data  is  then  graphically  fit  with  aid  of  a 
Houston  plotter  to  another  Pearson-IV  distribution.  From 
this  curve,  another  corrected  distribution  is  obtained. 

The  procedure  is  repeated  until  the  assumed  profile  is 
stable,  and  the  reverse  corrected  profile  agrees 
sufficiently  with  the  original  measured  data.  "Agrees 
sufficiently"  might  be  taken  to  mean  that  the  reverse 
correction  lies  within  the  experimental  error  of  the 
original  measured  data. 


V.  RESULTS  AND  DISCUSSION 


The  results  of  the  corrections  for  each  dose  is  shown 
on  graphically  in  subsequent  pages.  There  are  three  plots 
for  each  profile.  The  first  page  (Fig  (a))  shows  the 
measured  data,  plotted  as  circles,  along  with  its 
Pearson-IV  fit  .  Figure  (a)  also  shows  the  final  corrected 
Pearson-IV  distribution,  which  represents  the  real  profile. 
The  depletion  widths  calculated  from  this  profile  are  used 
to  correct  the  measured  data,  which  are  then  plotted  as 
triangles,  scattering  around  the  corrected  Pearson-IV 
profile.  The  depletion  widths  are  also  used  to  correct  the 
mobilities,  shown  in  the  bottom  of  the  graph  as  squares. 

The  second  page  (Fig  (b))  shows  the  original  measured 
data  and  the  reverse  corrected  profile  from  the  assumed 
real  profile.  In  each  of  the  doses,  the  reverse  corrected 
profile  sufficiently  fits  the  measured  data,  so  that  we  may 
conclude  that  the  iteration  procedure  has  indeed  produced 
the  real  profile,  and  thus,  the  correction  method  is  a 
success.  Convergence  of  the  assumed  real  profile  to  the 
final  true  profile  usually  required  approximately  7 
iterations  of  the  correcting  procedure.  However  it  was 
found  that  lower  doses  were  more  difficult  to  correct.  The 
reason  for  this  is  that  the  initial  depletion  width  takes 
up  more  of  a  substantial  part  of  the  profile.  Thus,  when 
the  newly  created  data,  which  extends  from  some  initial 


depletion  width,  is  fit  to  a  Pearson-IV  curve,  there  is 
more  uncertainty  of  how  well  the  Pearson-IV  fits  the  true 
profile  within  the  depletion  width.  There  is,  however,  no 
uncertainty  about  whether  or  not  an  assumed  Pearson 
distribution  is  the  actual  real  profile,  since  a  reverse 
calculation  will  immediately  reveal  that  it  either  does  or 
does  not  match  the  original  measured  data. 

On  the  final  page  (Fig  (c))  is  plotted  the  corrected 
data,  along  with  the  assumed  the  Pearson-IV  fit,  the  C-V 
data  and  finally  the  LSS  distribution  for  the  particular 
dose.  All  three  should  be  fairly  close  as  each  represents 
the  real  profile.  This  is,  however,  not  the  case.  In  each 
case  the  C-V  profile  is  below  the  corrected  Hall  profile, 
and  both  the  C-V  and  Hall  profiles  are  much  broader  than 
the  LSS  with  much  lower  peaks. 

By  integrating  the  corrected  profile  to  the  depletion 
width  we  can  find  the  approximate  density  of  surface 
states.  If  this  value  is  added  to  the  original  sheet 
concentration,  a  better  approximation  of  the  electrical 
activation  efficiency  may  be  obtained.  Table  I  summarizes 
the  number  of  surface  states,  the  previous  activation 
efficiency,  and  the  corrected  activation  efficiency  for  all 
the  samples  below.  From  the  apparent  and  real  activation 
efficiencies  it  is  obvious  that  the  higher  the  dose  the 
less  drastic  are  the  changes  in  the  measured  and  real 
profiles,  as  expected.  However,  note  that  for 


TABLE  I 


ACTIVATION  EFFICIENCY  AND  SURFACE  STATE  DENSITY 


The  apparent  activation  efficiency,  from  the  initial 
measured  sheet  carrier  concentration,  as  well  as  the 
corrected  activation  efficiency,  due  to  the  addition  of 
the  surface  state  density,  are  presented  here  for  each 
dose  that  was  measured. 


DOSE  cm“^ 

SURF.  STATES,  cm"^ 

APP.  ACTIV. 
EFF.  (%) 

REAL  ACTIV 
EFF .  ( % ) 

2  X  10^^ 

1.06  X  10^^ 

25.51 

78.54 

4  X  10^^ 

9.41  X  10^^ 

39.62 

64.23 

6  X  10^^ 

7.78  X  10^^ 

27.89 

40.86 

8  X  10^^ 

9.06  X  10^^ 

25.56 

36.89 

1  X  10^^ 

8.73  X  10^^ 

21.74 

30.48 

3  X  10^^ 

1.38  X  10^^ 

35.79 

40.40 
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2  X  10  cm  the  activation  efficiency  increases  three 
fold.  It  is  also  noted  from  Table  I  that  the  density  of 
surface  states  appears  to  be  independent  of  dose^  with 

statistical  fluctuations  around  approximately 

11  -2 
9  X  10  cm 

Since  there  seemed  to  be  a  good  deal  of  discrepancy 

between  the  C-V  and  Hall  profiles,  the  Hall  data  and  C-V 

data  from  a  sample  with  high  electrical  activation  are 

compared  with  a  non-electrical  means  of  profiling,  namely. 

Secondary  Ion  Mass  Spectrometry  (SIMS).  The  Hall  and  SIMS 
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data  for  a  sample  of  fluence  6  x  10  cm  were  taken 
from  Yeo  et  al.  (17),  and  a  C-V  profile  was  made  to  compare 
with  this  data.  The  data  is  presented  in  the  same  form  as 
above,  in  figures  19  (a)  through  (c),  except  that 
mobilities  are  excluded  and  the  SIMS  profile  is  included. 
The  results  shown  in  figure  19  (c)  are  very  interesting. 

The  C-V  profile  has  almost  the  same  shape  as  the  corrected 
Hall  profile,  and  yet  is  much  smaller  in  magnitude.  On  the 
other  hand,  the  SIMS  data  has  almost  the  same  peak 
concentration  as  the  corrected  Hall  profile,  however  it  is 
shifted  towards  the  surface.  Erroneous  measurements  of  the 
sputtering  rate  for  the  SIMS  data  and  the  etching  rate  for 
the  Hall  data  could  possibly  account  for  this  slight  shift 
and  drop  in  peak  concentration. 


U  ^  O  Measured  Hall  Data 


^  ^  A  A  Corrected  Hall  Data 


Q  □  □  □  Corrected  Mobilities 


« - *• - * - K  C-V  Profile 


B - B - a - n  SIMS  Profile 


LSS  Profile 


Figure  13.  Symbology  for  Measured  and 
Corrected  Data 


VI.  CONCLUSION 


The  real  and  apparent  carrier  profiles  have  been  shown 
to  be  related  by  the  differential  expression  (18)  and  so  if 
the  depletion  width  is  known  as  a  function  of  depth  the 
problem  is  immediately  solved.  However,  if  only  the 
apparent  profile  is  known,  there  is  no  direct  way  of 
finding  the  depletion  widths,  or  the  real  profile.  The 
method  undertaken  here  has  been  to  assume  a  reasonable 
distribution  as  a  first  guess,  and  to  calculate  from  it  a 
set  of  depletion  widths.  These  depletion  widths  are  used 
in  conjunction  with  the  measured  profile  in  equation  (17) 
to  find  an  updated  or  next  approximation  to  the  real 
profile.  It  has  been  found  that  by  iterating  this  process, 
the  real  profile  attains  a  stable  configuration,  and  no 
longer  changes.  This  means  that  the  real  profile  is 
consistent  with  the  measured  profile  through  equation  (17). 
Thus  we  have  found  the  real  profile  and  its  associated 
depletion  widths.  Correction  of  measured  profiles  has 
been  found  to  take  approximately  7  iterations.  This  extra 
effort  is  worth  while  since  it  takes  several  hours  to 
gather  a  Hall  profile. 

It  was  comforting  to  find  that  the  corrected  Hall 
profiles  were  much  closer  to  the  C-V  profiles  than  were  the 
original  measured  Hall  profiles,  thus  giving  an 
experimental  verification  of  the  correction  procedure.  It 


was  also  observed  that  the  C-V  profiles  were  always  much 
smaller  in  peak  concentration  than  the  corrected  Hall 
profiles.  This  is  possibly  due  to  conduction  through  the 
barrier  at  these  peak  concentrations.  Finally,  both  the 
Hall  and  C-V  profiles  were  found  to  have  peak 
concentrations  much  lower  and  standard  deviations  much 
larger  than  those  predicted  by  the  LSS  theory. 

There  are  several  advantages  of  the  differential  Hall 

method  over  the  C-V  method.  For  one,  the  Hall  method  works 

for  profiling  p-type  materials  while  C-V  profiling  only 

works  for  n-type.  Another  advantage  is  that  the  Hall 

method  also  measures  the  mobility  profiles.  In  addition, 

C-V  profiling  is  not  suitable  at  higher  doping  levels,  such 
12  -2 

as  6x10  cm  ,  since  avalanch  breakdown  will  occur. 

Most  importantly,  the  real  carrier  profiles  can  be  obtained 
right  up  to  the  surface  from  the  measured  Hall  profile 
using  the  correction  method  described  here.  However,  the 
C-V  profiles  can  be  obtained  only  beyond  the  original 
depletion  width.  These  advantages  make  differential  Hall 
profiling  an  important  alternative  to  the  simpler  and 
faster  C-V  profiling  and  necessitate  the  use  of  the  above 
procedure  to  correct  the  measured  Hall  profile. 


The  four  constants  a,  b^,  bj^f  b2f  of  the 
differential  equation  (3)  are  defined  in  terms  of  the  four 
moments  of  the  distribution: 

Rp  •J' Kf  (x)dx 


2 


./u-v 


(x)dx 


(x) dx 


^f(x)dx 


Where  R  is  the  projected  range,  ^  is  the  standard 
P  P 

deviation,  y  is  the  skewness,  and  P  is  the  kurtosis.  In 


terms  of  these,  the  four  constants  become: 


where 


A  =  10/9  -  12y^  -18. 

Equation  (3)  may  now  be  integrated  to  obtain  an  analitic 
function,  in  terms  of  these  constants.  This  function  is 
listed  in  line  180  of  the  program  in  appendix  B. 

The  effect  of  these  parameters  on  the  distribution  is 
demonstrated  in  the  following  graphs,  where  only  one 
parameter  is  varied  at  a  time.  The  default  values  are  a=l 
Nmax“80000'  (7p=0‘2,  y=0,  and  /S=10.  The  effect  of 

kurtosis  is  not  demonstrated  since  it  is  very  small. 
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DHCOia 

l  oi 


nc/'' 


DAVID  H.  ELSAESSER 
1  SEPT  1985 


THK-^FEOGRAH/^CORBECTS  A  KEASUREO  HALL  PROFILE  BT  FINDINO  ^ 
AN  APPROPRIATE  AND  SELF-CONSISTENT  SET  OF  DEPLETION  WIDTHS  . 


hf  ((  er-^erc:'^'. 


y 


u  or  c> 


■iz 


\ 


DIN 


X(100),T(100).XE(100),T1(100),NS(100>.NST(100).  NH<100).XN(100).XR(100).NR(1 
00>,XSF(100),NRF<100),NSn(100),NnP(100).nX(30),HOB(30>,HH(30> 

25  DIM  H(100),HF(100).NSC(500). 01(100). F(10>,CX<100>,CN(100t 
CLS) CLOSE) LT-L06( 10 ) l CORR>0) REVC0RR«0 
OPEN*COH1!4BOO,E,7,1’  AS  #1 
PRINT  •!,•))  • 

OEF  FNVBS(N)<i.6-.0259«L06(4.352E+17/N) 

OEP  FNU(N)<i3.  8051E+07»<  (VBS-.02S9)/N)‘^.5 
AA«0 

GOSUB  2000 

DBF  FNV(r)«FrX(7400«(r-CC)/rE> 

DEF  FNH(X)>FIX(5000»(X-AA)/XR) 

G05U8  1020) GOSUB  10000) GOSUB  7000 
OEF  FNA(KT. SR)-10*RT-12tSK*SK-18 


30 

40 

50 

53 

57 

60 

70 

B5 

90 

95 

120 


'DEF  CONSTANTS 
'A,B0,B1,B2 


130  OEF  FNB0(SD,KT. SK)«-SD«SDt(4*RT-3*SE*SK)/A 
140  DEF  FNB1(SD.I;T,SK)«-SK*SD*(KT>3)/A 
150  DEF  FNB2(KT. SK)«-(2*KT-3»SK»SK-6)/A 
160  DEF  FND(B0,Bl,B2)«4i80»B2-Bl*Bl 
170  DEF 

FNN(X)»(LOO(NO)*LOG(ABS(BO+B1»(X-XN  +  B1)+B2»(X-XH+B1>''2) )/2/B2-ATM(  (2tB2t(X- 
XH+B1)*B1)/DH)«(B1/B2«2*B1)/DH)/LT  'DEFINE  LOO  FORM  OF  PEAR-IV 
175  DH>SQR(0) 

IBO  DEF 

FMEM(X)«M0«(ABS<80  +  81#(X-XN*B1)*B2«<X-XN*BI)^2) >^(l/2/B2>*EXP(-ATN(  (2»Br*(X 
-Xn+El)+Bl)/DH)»(B1/B2*2»B1)/DH)  '  DEPJNE  PEAR-IV 

185  INPUT  'WANT  TO  INPUT  PREVIOUS  DATA)  '.ANSOIF  ANSO'T*  THEN  4000  ELSE 
GOSUB  7300 


200 

210 

230 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 


PRINT 

PRINT 

PRINT 

PRINT 

INPUT 

INPUT 

INPUT 


SD>.06)SX>.5iKT>10 
GOSUB  14000 
GOTO  260 

INPUT  "WANT  TO  TRT  ANOTHER  ITITAL  SET  OF  PARAMETERS)  '.ANSI 
IF  ANSf^’N'  THEM  440 
PRINT  ‘STANDARD  DEVIATION)  ■;SD 
‘KURTOSIS)  ‘jRT 
’SKEWNESS)  ‘iSK 
‘PEAK  POSITION)  ’)Xn 
‘MAX  IN  CONCENTRATION)  *>Nn 
'NEW  STAND.  DEV.)  ’.TSD)IF  TSDOO  THEN 
‘NEW  KURTOSIS)  ‘.TKT.'IF  TKTOO  THEN  KT" 

‘NEW  SKEWNESS)  '.T3R)IF  TSKOO  THEN  SR- 


SD-TSD 

■TKT 

■TSK 


71 


% 


3*0  INPUT  'NEW  PEAK  POSITIONi  •.TXn:IF  TXHOO  THEN  X(1-TXB 

3S0  INPUT  "NEH  MAX  IN  CONCENTRATION*  •,TNN»IF  THMOO  THEN  HH«TNN 

360  INPUT  ‘REPLOT  AXES  (NO-0,  TES-1)  i  *,C001: INPUT  ‘PRINT  PARANETEBS ( NO-0 

YES«1>  ‘.CODE 

370  60SUB  1*000 

370  IF  C001«0  THEN  *10 

*00  GOSUB  lOOOtGOSUB  7200 

*10  IF  C0&2-1  THEN  GOSUB  3000 

*15  CODl>0tC0O2>0 

*20  GOSUB  5000 

*30  GOTO  2*0 

*40  GOSUB  16020 

*50  SOO-SSi KT0>KTi SK0>SKi SI)B«S0i  KTB-KT:  SKB>SKi XHE-XN: HKB«NN 
*55  VBS>FNVBS(NH)iH-FNH(NH)iXH«Xn'»^H:FOR  K>1  TO  NE:  XR  ( K )  >XR  ( K)  4^U:  NEXT 
*65  INPUT  ‘ARE  YOU  FINISHED  COREECTING*  ‘.ANSJsIF  ANS»»‘Y’  THEN  20000 
*70  GOSUB  13*501 GOSUB  15000 
*75  GOTO  *65 
1000  ' 
loot  ' 

1002  ' 

1003  '  DRAW  AXES,  PUT  LABELS  ON 

100*  ' 

1005  ' 

1006  ' 

1010  GOSUB  2000 

1020  INPUT  ‘PRINT  AXES  HITH  PEN*  i  ‘.PEN 
1030  PRINT  *1,*  EC  1  A  1*20,2000  0  "I 
1035  PRINT  *1,‘  P'tiPRlNT  *1,  USING  ***tPEN 

10*0  PRINT  A  U  0,0  D  B  5000,0  0,7*00  -5000,0  0, -7*00  5000,0  0,7*00 

-5000,0  0,-7*00‘j 
1060  PRINT  •1,‘A  0,0  * 

11*0  ' 

11*1  ' 

11*2  '  OBAN  r  TICS  ON  CONCENTRATION  AXIS 

11*3  ' 

11*5  ' 

11*6  ' 

1150  Cl»FIX(CCC)iDl»FIX(DDC*.l)iN»0 
1160  FOR  J’«C1  TO  D1 
1170  FOB  K»1  TO  9 
1180  Y«LOG(K)/LT+JsL«FNV(Y) 

1190  IF  Y<(CCC-.05)  THEN  1250 
1200  IF  YXDDC  +  .05)  THEN  1300 
1210  IF  K-1  THEN  12*0 

1220  PRINT  #1,‘A  ■jHj‘,‘:L»‘D  R  75.0  -75.0  U*| 

1230  GOTO  1250 

12*0  PRINT  #1,‘A  ‘inj*,‘lL»‘D  B  150,0  -150,0  U  -500,-100  S13* ‘ : STR* ( 10 > ;  *. 
‘j*  -50,100  ‘j‘S13+‘>STH»(J)j‘  ‘I 
1250  NEXT  K 
1260  NEXT  7 

1300  '  NON  ORAN  TICS  OH  OTHEB  SIDE 
1310  YR-YRn>OD«ODH:CC>CCN 

1350  Cni«FIX(CCn+. 05)sDni«FIX(DDW+.  1)iN«5000 
1360  FOR  J«Cni  TO  OHl 
1365  PRINT  J 
1370  FOR  K«l  TO  9 


1375  PRINT  K 

1300  Y»LOG(K)/LT*J:L»FNV(Y) 

1390  IF  Y<(CC«-.05>  THEM  1450 
1400  IF  Y>(DDn+.05>  THEM  1470 
1410  IF  K-1  THEN  1440 

1420  PRINT  *1,’A  R  -75,0  75,0  U*s 

1430  GOTO  1450 

1440  PRINT  #1,‘A  ‘jMj’.'fLj’D  R  -150.0  150,0  "is  IF  COD3»0  THEN  1450  ELSE 
PRINT  •!,•  U  50,-100  S13+ * j STRl ( 10 > j *  •»*  -50,100  * » -SIR* * ; STR» ( J ) ;  *  • 

1450  NEXT  K 
1460  NEXT  7 

1470  YR«YRCiCC>CCC:00>0DC 
1480  PRINT  il,  ■  U* 

1500  ' 

1501  ' 

1502  ' 

1503  '  ORAM  X  TICS  ON  DEPTH  A^IS 

1504  ' 

1505  '  FIRST  ON  THE  BOTTOM 

1506  ' 

1520  L«Oi TEST-0 

1530  FOR  X-.05  TO  BB«.025  STEP  .05 
1540  H-FMH(X) 

1550  IF  TEST-1  THEN  1580  'GOTO  BIG  TIC  HARK 

1560  PRINT  41, -A  R  0,50  0,-50  U*;«  TEST-1  'SHALL  TIC  HARK 

1570  GOTO  1590 

1500  X«F1X<10»X)/10IPRINT  #1, 'A  ‘i  Uf  %  ‘flf  "J>  R  0,  100  0,-100  U  -150,-200 
S13>*( ST8$(XI) •.  ’iiTEST-O 
1590  NEXT 

1595  PRINT  #1,  *A  10000,10000  •» 

1720  L-74001 TEST-0  'NON  PLOT  TICS  ON  TOP 

1730  FOR  X-.05  TO  B8«.025  STEP  .05 
1740  H-FNH(X) 

1750  IF  TEST-1  THEN  1780  'GOTO  BIG  TIC  HARK 

1760  PRINT  #1, "A  •}H»«i'iL»‘0  B  0,-50  0,50  U-jtTEST-l  'SHALL  TIC  HARK 
1770  GOTO  1790 

1780  PRINT  #1, -A  ’jH? •, *(11  "D  B  0,-100  0,100  G  *»» TEST-0 
1790  NEXT 

1795  PRINT  #1,  -A  10000,10000  •» 

1800  RETURN 
2000  ' 

2001  ' 

2002  ' 

2003  '  INPUT  THE  PLOTTER  PARAMETERS 

2004  ' 

2005  ' 

2010  INPUT  'MAX  IN  OEPTHi  '.TSBUF  TEBOO  THEN  BB-TBB 

2020  INPUT  -MIN  IN  COHCENTHATIONJ  •,TC«IF  TCOO  THEN  CCC-LOG ( TC ) /LT 

2030  INPUT  -HAX  IN  CONCENTRATION*  -.TDslF  TDOO  THEM  DDC-LOO (TD I /LT 

2035  INPUT  -MIN  IN  CONCENTRATION  PLOTTING*  ■,THP*  IF  THPOO  THEM 
up-log (TMP I /LT 

2040  INPUT  -HIM  IN  MOBILITY*  ‘.TCCHiIF  TCCHOO  THEN  CCH-LOGCTCCH) /LT 

2050  INPUT  "MAX  IN  MOBILITY*  -.TODHSIF  TDDHOO  THEN  DDH-LOGI TDDH I /LT 

2060  INPUT  -UANT  TO  PUT  MOBILITY  SCALE  ON  (1-YES.  0-N0)i  -.CODS 
2080  XS-(BB-AA)i YRC- ( DDC-CCC ) : YRH- ( DDH-CCH ) *  YE-YRC* CC-CCC*  DO-DDC 
2090  IF  ABSCYRC-YHHO.OOl  THEN  2010 


2100  SETUBN 

3000  ' 

3001  ' 

3002  ' 

3003  ' 

3004  '  PRINT  THE  PEARSON  IV  PASAHETESS 

3005  ' 

3006  ' 

3007  ' 

3008  ' 

3009  ' 

3010  INPUT  -PRINT  PARAMETERS  HITH  PEN*  t  *,  PEN: PRINT  *1.-  p-;t PRINT  #1, 
USING  -•‘;PEN 

3020  PRINT  *1,  -A  8000,5600  S14  X  MAX-  * 5 STRt ( XM ) ; ■ _  -j 

3030  PRINT  *1,  -A  8000,5400  S14  N  MAX-  *tSTR*(NH>t‘  ‘t 

3040  PRINT  *1,  -A  8000,5200  S14  STAND.  DEV.-  * ( STR$ ( SD ) ; ■_  -f 

3050  PRINT  *1,  -A  8000.5000  S14  SKEWNESS-  * ; STRt ( SK ) » ■_  “j 

3060  PRINT  #1,  -A  8000,4800  S14  KURTOSIS-  * > STR» ( KT ) ( *.  -j 

3080  RETURN 

4000  ' 

4001  * 

4002  ' 

4003  '  INPUT  EXISTING  SATA 

4010  FILN2*-PIL:it  +  *.DHC-:FII.N3i*FII.M**‘.COP-i2-0:  COP-1 

4020  OPEN  "I-,3, FILN2t:OPEN  -O’.I.FILNSt 

4030  INPUT  *3,  NE,XH, NR, SOO.SRO.KTO 

4040  FOR  J-1  TO  NE: INPUT  *3, XB(T) , NR< J ) t NEXT  J 

4050  IF  E0P(3)  THEN  4100 

4055  PRINT  XH,NM,SDO,SKO.KTO: INPUT  “GO  FROM  HERE:  •,AMS»:IF  ANS»«*Y*  THEM 
4100 

4060  PRINT  •2,NE,XM.Mn. SDO.SKO.KTO 

4070  FOR  J»l  TO  NE:PRINT  *2, XH( J) , NR (J) : NEXT 

4080  IF  Z-0  THEN  XMB-XRt NRB-NR: SOB-SSO:  SKE-SKO:  KTE-KTO 

4090  Z«Z*l:GOTO  4030 

4100  CLOSE  3: KILL  FILNZt 

4120  CORR-Z 

4130  INPUT  -ARE  YOU  FINISHED  CORRECTING:  *, ANSI: IF  AHSI-’Y*  THEN  4140 
4135  IF  Z-0  THEN  240  ELSE  468 

4140  PRINT  #2. NE,XH, NR. SDO.SKO.KTO: FOB  T-l  TO  NE:PBINT  42. XB( J) . NR ( I > :  NEXT 
4150  GOTO  20000 
5000  REH 

5005  INPUT  -PLOT  PEAH4  HITH  PEN*  :  •,PEM:PBINT  41,*  P-(:PRIHT  41,  USING 
*4-jPEM 

5010  Y»FNN(AA):M»FHH(AA1:L-FNV(Y):  IF  YCMP  THEN  5030 
5020  PRINT  41, -U  A  •:M:*,-»L|*  D  *: 

5030  DX-XR/100 

5035  IF  PEN-0  THEN  RETURN 

5040  FOR  X-AA  TO  BB  STEF  DX 

5050  H»FNH(X): Y-FNN(X):L»FNV( Y): IF  Y<»MP  THEM  PRINT  41,  ■  U  ’::GOTO  5060 
5055  PRINT  41. Mj"  -iLi*  D’|: 

5060  NEXT 

5070  PRINT  41, -U  10000,10000  *» 

5100  RETURN 

7000  REM 

7001  RER  CALCULATE  THE  PROFILE 

7002  REH 


7003  BEN 

7010  INPUT  ‘INPUT  THE  FILENAHE  CONTAINING  HALL  PROFILE  DATAi 
•,FILHt:FILNl»»FILN*+‘.DAT* 

7030  OPEN  "I'.B.FILNl* 

70AO  INPUT  #2, I, HAT 
704S  TT«0«XE(0)"0 
7050  FOB  K-1  TO  I 
7060  INPUT 

*2. TI (K), NN(IC): TTaTT^TI (K ) s  XE <n ) «TT«BAT/10000:  DI ( K ) «TI ( K I cBAT/lOOOO 

70B0  NEXT 

7085  CLOSE  *2 

7087  NN-0 

7090  FOB  K»l  TO  I 

7100  XH(K)«<XE(K)+XE(K-in/2 

7105  NB(K)>NH(K)iXB(i:)«Xn(Kls IF  NB(R>>NH  THEN  NH>NB ( K ) > XH>XB ( K) 

7110  NEXT  K 
7120  NE»I 

7130  OPEN  •1*,2,F1LNI+*.N0B* 

7135  INPUT  *2,  IN 
7140  TOTN«0:NX<0)-0 
7145  FOB  K-1  TO  IN 

7150  INPUT  #2,X,NOB<K)JTOTN«TOTN*XsNX(K»«TOTH»BAT/10000 

7153  NEXT 

7160  CLOSE  *2 

7200  HEN  PLOT  THIS  DATA 

7202  INPUT  ‘PLOT  DAT,^  HITH  PEN#  «  *,  PEN:  PRINT  #1,  •  P*»t  PRINT  *1,  USING 
‘i'jPENj 

7210  FOB  K«1  TO  I 

7215  IF  NN(K)<0  THEN  7240 

7220  X«FNH(XN(K))tT*LOG<HH(Kl)/LT!lF  Y<HP  THEN  7240  ELSE  Y»FNV(y) 

7223  IF  X05000  THEN  7240 

7230  PRINT  #1,  ‘CC  ‘ ( X( * , * j Yj •33* j 

7240  NEXT 

7250  PRINT  #1,  •  U  A  10000. 10000‘j 
7260  RETURN 
7300  REN 

7310  FILN2t»FILN»*‘. DHC‘ 

7320  OPEN  ■0*,2,FILN2» 

7400  RETURN 

7500  HEN  PLOT  THIS  DATA  USING  TBIAirLES 

7502  INPUT  ‘PLOT  DATA  HITH  PEN#  «  •,PEN:PBINT  #1.*  P‘i: PRINT  *1.  USING 
"••;PEN« 

7510  FOB  K»1  TO  NE 
7513  IF  NR(K)<0  THEN  7540 

7320  X«FNH(XH(K) ): Y»LOO(NR(K) )/LTt IF  Y<HP  THEM  7540  ELSE  Y»FNV(Y) 

7523  IF  X03000  THEN  7340 
7530  PRINT  #1,  X(‘.*;Yj‘  H44  ‘i 
7540  NEXT 

7530  PRINT  #1,  •  U  A  10000, lOOOO'i 
7360  RETURN 
7600  REN 

7610  FILH2»»FILNS+*.DHC‘ 

7620  OPEN  ‘O’.R-.  FILN2$ 

7700  RETURN 


8002  '  CALCULATE  THE  COERECTEO  PROFILE 

8003  '  FINS  DEPLETION  UIDTHS 

8004  ' 

8003  ' 

8010  '  COMPUTE  SHEET  CONCENTRATIONS 
8020  PRINT  ■INTEGHATINO  N<X>' 

8040  INPUT  "XMAX  FOR  INTEGRATION!  ‘.XMAX 

8030  HO>XHAX/SOO 

8033  NSC(0)>0 

8060  FOR  J-1  TO  300 

8070  F1«FNEN( ( J-. 3» *H0 > ! NSC < J ) «NSC ( T-l )+Fl*H0/10000 
8073  PRINT  T.NSCCJ) 

8090  NEXT 

8093  H1»(XHAX-.1)/100 
8100  INPUT  *U0i  *,H(0) 

8110  FOR  J>0  TO  100 
8113  Z-0 

8120  IF  JOO  THEN  U(J')-U(J-1) 

8130  X«(H( J)*J»Hl)/HOsK«FIX(X)iDK»X-Kt Y«J*Hl/H0iL»riX(Y):DL«T-L 
8140  IF  X>S00  THEN  NE«J-l!GOTO  8230 

8130  N«INSC(K)+(NSC(K4l)-NSC(K) )*DK-NSC(L>-(NSC(L«1>-NSCIL) ) (DL) tlOCOO/U ( 7) 
8160  VBS>FNVBS(N)!UF«FNU(N) 

8163  2«2+l!lF  2>30  THEN  INPUT  'UFi  •,HF!Z«1 

8170  IF  ABS(H<7)-UF)<.0001  THEN  8210 

8180  H(7)>(U( J)'fHF)/2 

8200  GOTO  8130 

8210  PRINT 

8220  NEXT 

8223  NE*100 

8230  IF  REVCORR*!  THEN  RETURN 

8233  GOSUB  9000 

8300  REN  CORRECT  THIS  PROFILE 

8310  TNH<0 

8313  PRINT  NE 

8320  FOR  7«1  TO  NE 

8340  X8( J)»( (2»J-l)tHl+H< J)+U( J-1) 1/2 
8330  N8( J)«NNP(J)/(1+(U( J)-U( J-1) )/Hl) 

8360  IF  NR(J)?TNN  THEN  TNM-NE ( J)  i  TXN-X8( J) 

8363  PRINT  J, XR ( J ) . NR ( J) 

8370  NEXT 

8300  REM  PLOT  THE  CORRECTED  CURVE 

8510  INPUT  "PLOT  DATA  HITH  PEN#  i  •,PEM! PRINT  #1,*  P*) SPRINT  #1,  USING 
■f;PEN» 

8513  IF  NRdXO  THEN  N8(i)-1E+13 

8520  X-FNH(XR(1) )s Y*L0G(NR<1) )/LTs IF  T<nP  THEN  8530  ELSE  Y-FNVIY) 

8323  PRINT  #1,  *U  •jXj'.'jYi’D  "j 

8530  FOR  IC-2  TO  NE 

8540  IF  NRCKXO  THEN  8570 

8330  X>>FNH(XR()C)  ) i  Y>FNV ( LOG(N8 ( 1C) ) /LT> 

8355  IF  X«>3000  THEN  8370 
8360  PRINT  #1.  Xj  *,  'i  Yi 

8563  IF  Y<«0  THEM  PRINT  #1,  •  U  *>  ELSE  PRINT  #1,  ■  D  •» 

8370  NEXT 

8380  PRINT  #1,  *  U  A  10000,10000  'j 
8590  RETURN 
9000  ' 


9020  FOH  J-l  TO  NE 

9030  X»( )/2iMnP<J)«FMEM<X»*ll»(U(J)-U(J-l) )/Hl> 

9040  NEXT 

9030  Xn*XHBiNH-NnB:5D«SDBiSR«SXBiRT>rvTS<60SUB  14000 
9070  FOB  7-1  TO  NE 

90S0  NHP(7)«IFNEN(  (7-. 8)  tHl ) 'fNNP (7)  > /2 
9090  NEXT 
9100  BETUBN 

10000  BEN  INPUT  C-V  FILE 
10080  INPUT  *DOSE:  ',OOSE 

10060  INPUT  'LSS  RPt  '.BPi INPUT  ‘LSS  STAN.  DEV. l  '.SIG 

10070  N6>DOSE*10000/SIO/2.3066283tiDEF  FNLSSt X)«NO*EXP ( - ( ( X-RP ) /SIG ) '2/2 ) 
10073  INPUT  ‘BEGINNING  HEASUREO  SHEET  CONCENTRATION!  ‘.NSHO 
lOOBO  PRINT  ‘HAX  IN  LSSi  *,N6 

10100  INPUT  ‘NAME  OF  C-V  FILE!  *,  CFILN*t CFILN»»CFILN«  +  ‘ . DAT* 

10110  OPEN  ‘I*,4,CFILNt 
10120  INPUT  M.CN 

10130  FOR  7-1  TO  CN: INPUT  *4,  CX(7) ,  CN (7) > NEXT 
10140  CLOSE  4 
10200  RETURN 
10280  ' 

10233  INPUT  ‘PLOT  C-V  DATA  WITH  PEN#  I  ‘.PENtPRINT  •!, ‘  P’pPRINT  tl,  USING 
‘♦‘;PEN| 

10260  L»FNH(CX(1I >i T-LaG(CH(ll )/LTin>FNV(TM IF  TCNP  THEN  10270  ELSE  PRINT 
*1,‘U  ‘jL»‘,‘jN;‘D  ‘j 
10270  FOR  7-1  TO  CN 

10280  L«FNH<CX(7))iT>L00(CN(7n/LTtK»FNV(TM  IF  TCHP  THEN  10300 
10283  IF  X73000  THEN  10300 
10290  PRINT  Ol.Li ‘, ‘jHj  ■,  1141  ‘ 

10300  NEXT 

10310  PRINT  *1, *U  10000,10000* 

10400  RETURN 

10300  '  PLOTT  LSS  FUNCTION 

10320  OX-XR/100 

10330  INPUT  ‘PLOT  LSS  HITH  PEN#!  *, FEN: PRINT  tl, ‘  P‘|i PRINT  »1,  USING 
‘•‘>PEN| 

10340  L-FNH(OMn-FNV(LOO(FNLSS(On/LT)!PRINT  *1,  ‘  U  *jLj‘,‘|Hf'D  •j!2«0 
10330  FOR  X-AA  TO  BB  STEP  DX 
10360  Z-Z+1 

10370  L>FNH(X)i T-LOG(FHLSS(X) )/LTi  n«FNV(TI 

10373  IF  T<-nP  THEN  PRINT  *1,  ‘U  *! X-BB+l! GOTO  10600 

10380  IP  Z-4  THEN  PRINT  #1, ‘  U  ‘)L|‘,‘?Mj‘  D‘;!Z»OiGOTO  10600 

10390  PRINT  #1,  L| *, ‘»Bj 

10600  NEXT 

10630  PRINT  01,*  U  10000.10000*1 
10660  RETURN 

lOSOO  '  PLOTT  THE  CORRECTED  NOBILITIES 

10803  INPUT  ‘PLOT  NOBILITY  HITH  PEN#  l  ‘.PENsPRIHT  *1,*  P‘t!PHIHT  #1,  USING 
‘•‘sPENj 

10810  CC-CCHiDO-OOHiYR-rRN 
10820  FOR  7-1  TO  NEH 

10830  X-(HX(7) tHX(7-l)+HH(7)»HH(7-l»  >/2!X»FHH(X> 

10833  IF  HOB(7)<-0  THEN  10870 
10840  Y-L00(H0B(7) )/LT 
10830  IF  T<CCH  OR  T>CCn+l  THEN  10870 
10833  IF  X:3000  THEN  10870 


10860  T»FNV(Tl:PHINT  #1,  1142*| 

10870  NEXT 

10880  CC«CCC:D0«0DC: YS-TRC 
lOYOO  BETUSN 
11000  ‘ 

11100  INPUT  “NAME  OF  DATA  FILEi  ■,  CFH-H»« CFXLN»»CFILH»* ■ . DAT" 

11110  OPEN  ■I',4,CFtLN« 

11120  INPUT  *4, ON 

11130  FOE  J-1  TO  CN:INPUT  #4, X( J) , Y < J) « NEXT 
11140  CLOSE  4 

11235  INPUT  "PLOT  PBOFILE  WITH  PEN#  »  ‘.PENsPRINT  *1,*  P'jsPBINT  *1,  USING 
■••(PEM| 

11260  L-FNH(X(l))«Y«LOG<Y(l))/LTiN-FNV(T)»IF  T<BP  THEM  11270  ELSE  PBINT 
#1,*U  •»Lj*,*jN5"D  •» 

11270  FOB  T-1  TO  ON 

11280  L«FNH<X<J)UY»LOQ<T(J))/LTsn-FMV(Y)iIF  YCNP  THEN  11300 
11283  IF  L03000  THEN  11300 
11290  PBINT  *l,Lj  "i  *f  M»  *.  W42  • 

11300  NEXT 

11310  PBINT  «1,*U  10000,10000* 

11400  BETUBN 
13430  BEN 

13480  PBINT  ’STANDARD  DEVIATION*  ’iSOO 
13490  PRINT  'KUBTOSISi  ’iKTO 
13300  PRINT  'SKEUNESSt  •;SR0 
13310  PBINT  'PEAK  POSITION*  "iXN 
13520  PRINT  *KAX  IN  CONCENTRATION*  **NN 

13330  INPUT  "NEH  STAND.  DEV.*  '.TSEJlF  TSDOO  THEN  S00»TS0 

13340  INPUT  *NEH  KUBTOSIS*  '.TKTiIF  TKTOO  THEN  KTO-TKT 

13530  INPUT  'NEH  SKEHNESS*  ’.TSKJlF  TSKOO  THEN  SK0»TSK 

13360  INPUT  'NEH  PEAK  POSITION*  ".TEN*  IF  TXNOO  THEN  Xn»TXB 

13370  INPUT  'NEH  MAX  IN  CONCENTRATION*  *,TNM»IF  TNNOO  THEN  MB»TNH 

13380  INPUT  "BEPLOT  AXES  (O-NO,  YES»1>  »  ‘.CODl* INPUT  "PRINT  PABABETEBS 

<N0>0,  YES-1)  *  ’,COD2 

13390  INPUT  "CLEAB*  *,01* IF  CL-1  THEM  PBINT  #1, *Z  »**» 

13600  IF  CODl-O  THEN  13620 

13610  GOSUB  lOOOiGOSUB  8300 

13620  SD-SOOsKT-KTOsSK-SKO* GOSUB  14000 

i:.623  IF  C002-1  THEN  GOSUB  3000:COD2-0 

13630  SDO«SO*KTO«KT:Sr.O«SK 

13640  BETUBN 

14000  BEN 

14020 

A-FMA(KT,SK)*B0-FHB0(SD,KT,SK)*B1-FNB1(SD,KT, SK>iB2-FKB2(CT,SK)i0»FND(B0,Bl 

B2 ) 

14080  IF  D<-0  THEN  PRINT  "OISCBIBINANT  ■  *10*00X0  14100 
14090  GOTO  14250 
14100  BEN 

14110  PRINT  ‘INPUT  HEN  PABABETEBS  ' 

14120  PRINT  ’STANDARD  DEVIATION*  '{SO 
14130  PRINT  "KURTOSIS*  * !  KT 
14140  PRINT  'SKEHNESS*  *tSK 
14130  PRINT  "PEAK  POSITION*  "JXB 
14160  PRINT  'BAX  IN  CONCENTRATION*  "*NH 

14170  INPUT  "NEH  STAND.  DEV.*  '.TSDslF  TSDOO  THEN  SD-TSD 
14180  INPUT  "NEH  KURTOSIS*  *,TKT*IF  TKTOO  THEN  KT-TKT 


14190  INPUT  ‘NEH  SKEHNESSt  '.TSKilF  TSKOO  THEN  SN«TSK 
14200  INPUT  ‘NEH  PEAK  POSITIONi  •,TXH:IF  TXNOO  THEM  XN-TXH 
14210  INPUT  'NEW  MAX  IN  CONCENTSATIONt  *,TNn:  IF  TNHOO  THEN  NN«TNII 
14220  GOTO  14020 

14250  0H>SaR<0)  I  N0>li M0aNH/FNEN(XH) 

14260  RETURN 
15000  ' 

15100  ' 

15110  PRINT  ’CHANGE  PAPER  'sGOSUB  5000 

15120  INPUT  *S0  YOU  HANT  TO  CORRECT  NOUi  *,ANS* 

15130  IF  ANSt>*T*  THEN  GOSUB  16000  ELSE  RETURN 
15132  TTXN-XNiTTMN-NN 
15135  GOSUB  eoOO 

15137  INPUT  *HANT  TO  PLOT  REVCORRECTEfi  OATAt  ’.ANStilF  ANSt«’T’  THEN  GOSUB 
17000 

15138  XH-TXN:NH-TNN 

15140  GOSUB  1000; GOSUB  8500; GOSUB  10250 
15150  INPUT  ’CHANGE  PAPER  THEN  RETURNi  ’.ANS* 

15160  INPUT  ’HANT  TO  PLOT  REVCORRECTEO  OATA>  ’.ANSI; IF  ANSt>'T’  THEN  GOSUB 
17000 

15300  RETURR 
16000  REN 

16010  SO-SOOs  SK>SK0i  KT«RT0t  C0BR>C0RR'»1 

16020  LFRINT  ’CORRECTION  *’>0088 

16030  LPBINT  ’XHAXt  *tXH 

16040  LPBINT  ’NHAXI  ’{NH 

16050  LPBINT  ’STANDARD  DEVIATION!  ’{SO 

16060  LPBINT  ’SREUNESSt  ’{SR 

16070  LPBINT  ’KURTOSISi  ’ | KTt LPBINT! LPBINT 

16080  PRINT  t2,NE,XM,NN,SD,SR,KT 

16090  FOR  J«1  TO  NEiPRINT  *2,  XR(T),  NR(T)t NEXT 

16100  RETURN 

17000  ' 

17100  '  CALCULATE  THE  REVERSE  CORRECTED  NEASUREO  DATA 
17110  INPUT  ’CHANGE  PAPER,  THEN  RETURNi ’, ANS» 

17120  GOSUB  1000 

17140  PRINT  ’NON  PLOTTING  ORIGINAL  HEASURED  DATA 
’:XH>XHB!Nn«NHB!SD-SDBiSK«SRBiRT-RTBiG0SU8  1400: GOSUB  5000 
17145  Nn>TTNHiXH-TTXn:SD«SDOiKT<KTOtSK-SROi GOSUB  14000; FOR  J*l  TO 
NE!X(J)>XR(J)iT(T)»NR(T)!NEXT 
17150  FOB  J»1  TO 

NE!X»( (2tJ-l)«Hl»(H(J'l+H(J-l) ) )/2iNR(J»»FNEM(X)#(l*(H(T)-H( J-1) >/H1)!XRIJ)» 
(0r-.5)*Hli  NEXT 

17170  PRINT  ’NOH  PLOTTING  REVERSE  CORRECTED  DATA  ’ 

17180  GOSUB  8500 

17190  FOB  J«1  TO  NE:XR(J')»X(a’)tHB<J’)»Y(J»iMEXT 
17300  RETURN 
20000  ' 

20001  ' 

20002  ' 

20003  ' 

20004  ' 

20005  '  FINISHED  CORRECTING  —  NOW  CHECK  FOR  SELF  CONSISTENCY 

20006  '  AND  PRINT  FINAL  PLOTS 

20007  ' 

20008  ' 


20009  ' 

20010  NEL-NE: XnLaXHi NHL-HHi  SDL«S00:  SI:L«SR0: KTL>KT0 
2001S  FOR  J*l  TO  NEL:XRF(J)*XR(J):NRP( J>*NR(J): NEXT 
20020  SS>SOO:SK*SXOiKT-KTO 

20030  GOSUB  14000 

20040  CLOSE  2iXF  COF-1  THEN  NAME  FILN3*  AS  FILH24 
200S0  PRINT  ‘INTEGRATING  FOR  FINAL  TIME*  t REVCORB>li GOSUB  8040 
20100  '  CALCULATE  THE  REVERSE  CORRECTED  MEASURED  DATA 
20110  INPUT  'CHANGE  PAPER.  THEN  RETURN! *, ANS« 

20120  GOSUB  1000 

20130  FOR  T-l  TO  H XR ( J) «XH ( J ) !  NR ( J > >NM(3) i NEXT! ME«I 
20140  PRINT  'NON  PLOTTING  ORIGINAL  MEASURED  DATA  ‘i GOSUB  7200 
20143  NE-NEL:XH>XHLiNH«NHL!  SO>SDL:SK>SKLiKT*RTLt GOSUB  14000 
20130  FOR  J-1  TO 

NE!X-< (2•J-l)•Hl♦(H(T)«U(;-ln)/2tNR(T)«FN£N(X)•(l♦(H(J)•U(J•l))/Hl)lXR(3)■ 
(J-.3)»HliNEXT 

20170  PRINT  ‘NON  PLOTTING  REVERSE  CORRECTED  DATA  * 

20180  GOSUB  8300 

20300  '  CALCULATING  DEPLETION  WIDTHS  FOR  MEASURED  DATA 

20303  PRINT*  NON  CALCULATING  DEPLETION  WIDTHS  FOR  MEASURED  DATA* 

20310  INPUT  *H0t  ‘,UF(0> 

20320  FOR  J>0  TO  I 
20340  Z-O 

20330  IF  TOO  THEN  UF( J)>HF(7>1> 

20360  X»(HF<JJ+XE(JI)/HOtK«FIX<X»iDK«X-K!r»XEtJ)/HO!L«FIX«T»!DL«r-l 

20370  IF  X>300  THEN  NE-J-HGOTO  20300 

20380 

N«<NSC(K)«(NSC(K«1)-NSC(R) ) •DR-NSC (L) - (NSC(L*1 ) >NSC (L> ) tDL) tlOOOO/WF ( 

20383  NSH( J)>NSC<IO*(NSC(C«l)-NSC(K)>*DK 
20390  VBS»FNVBS<N)IHT«FMH(N» 

20400  2«2*liIF  Z>30  THEN  INPUT  *HFi  ’.WTiZ-l 
20410  IF  ABS(HF(J‘)*UT)<.0001  THEN  20430 
20430  HF(J’)«(WF( J)>UT)/2 
20440  GOTO  20360 
20430  PRINT  7, HF(J) 

20460  NEXT 
20463  NE-I 
20470  GOSUB  23000 
20300  ' 

20310  X»HP ( P ) /HO! K»FIX( X) :  DK-X-K 

20340  LPRINT!LPRINT  ‘ORIGINAL  DEPLETION  WIDTH!  ‘,HF(0) 

20330  LPRINT  ‘APPARENT  ACTIVATION  EFFICIEHCTi  *( NSMOtlOO/DOSE 
20560  LPRINT  ‘TRUE  ACTIVATION  EFFICENCTi 
‘I (NSMO«NSC<K)«^(NSC(Rtl)-NSC(R) )tDR)*100/D0SE 
20570  LPRINT  ‘THE  APPROXIHATE  NUMBER  OF  SURFACE  STATES  IS! 

‘f  NSC(K)«’(NSC(R*1)>NSC(R)  I  tDR 

20600  '  PLOT  THE  MEASURED  DATA  AMO  CORRECT  THE  MEASURED  DATA 
20610  PRINT  ‘READY  TO  PLOT  MEASURED  DATA'! INPUT  ‘CHANGE  PAPER,  THEN 
RETURN', ANS$ 

20620  GOSUB  1000 

20630  FOR  7-1  TO  l! XR<7) ■XH(7) !  NR <71 -HN(7 1 ! NEXT 
20640  GOSUB  7200 
20630  FOR  7-1  TO  I 

20660  NR<7)-NH(7)/(1+(HF(7)-HF(7-1) >/DI<7) )!XR(7)-XM(7)«’(HF(71*HFt7-l» )/2 
20670  NEXT 

20675  PRINT  "PLOT  THE  CORRECTED  DATA  ‘ 


80 


20680  GOSUB  7500  'PLOT  TBI ANGLES  FOB  COBBECTEO  DATA 

20700  '  NON  PLOT  THE  HEASUBEO  AND  COB8ECTED  PEABSON  IV 

20710  PBINT  "PLOT  THE  PEABSON  IV  FIT  TO  THE  OBI6INAL  DATA* 

20720  Xn-XnB<NH-NHBiSD-SDCtSrv-SKBsfCT>KTB:  GOSUB  lAOOOi  GOSUB  5000 
20730  PBINT  ‘PLOT  THE  PEABSON  IV  FIT  TO  THE  COBBECTED  DATA* 

20740  XH-XHL:NN«NHLiSD*SDLiSR>SRLit:T«i:TLi  GOSUB  14000t  GOSUB  5000 
20742  INPUT  *UANT  TO  PLOT  MOBILITIES  *,AHS»ilF  ANS<«*T*  THEN  GOSUB  10800 
20745  GOSUB  lOOOtPBINT  ‘COBBECTEO  DATA  **60SU8  7500:PB1NT  ‘COBBECTEO 
FIT* s GOSUB  5000 

20750  PBINT  ‘NOU  PLOT  LSS  *: GOSUB  10500 
20760  PBINT  ‘NOU  PLOT  C-V  OATA‘l  GOSUB  10250 

20850  INPUT  ‘DO  YOU  UAHT  PLOT  AHOTHEB  HEASUBED  PBOFILEi  ‘.ANSI 
20860  IF  ANS4>‘T‘  THEN  GOSUB  11000 

21770  INPUT  ‘WANT  TO  PLOT  NOBILITIES  ‘.ANSdlF  ANS«-‘T‘  THEN  GOSUB  10800 
21850  INPUT  ‘DO  YOU  WANT  HOBS  COPIESi  *,ANS« 

21860  IF  ANS«<>‘N‘  THEN  GOTO  20110 
21904  END 

25000  '  CALCULATE  DEPLETION  WIDTHS  FOB  THE  NOBILITY  DATA 

25020  PBINT: FEINT  ‘CALCULATING  DEPLETION  WIDTHS  FOB  THE  NOBILITY  DATA 

*: PBINT 

25310  INPUT  ‘UOI  ‘,HM(0l 
25320  FOB  X-0  TO  IN 
25340  Z>0 

25350  IF  700  THEN  HN(7)«UH(7-1) 

25360  X>(HH(7)'^NX(7n/HO:K*FIX(X>iOR>X>KiT>HX(7)/HOtL>FIX(T):DL>T-L 
25370  IF  X>500  THEN  NEH-7-ltOOTO  25500. 

25380 

N« ( NSC <K)4<NSC(K>1) -NSC (K) )«OK-NSC(L)-<MSC <L>1) -NSC (L>  )#DL> *10000/85(7) 
25390  VBS>FNV8S(N)i HT«FNU(N) 

25400  Z«Z+l:IF  Z>50  THEM  INPUT  ‘HN«  ‘.WTtZ-l 

25410  IF  ABS(HH(7)-HT)<.0001  THEN  25450 

25430  UN<7)«(UH(7)*UT)/2 

25440  GOTO  25360 

25450  PBINT  7,HF(7I 

25460  NEXT 

25465  NEN-IN 

25470  BETUBN 
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